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The probability of simultaneous occurence of at least k spanning clusters has been studied by Monte Carlo 
simulations on the 2D square lattice at the bond percolation threshold pc ~ 1/2. The calculated probabilities for 
free boundary conditions and periodic boundary conditions are in a very good coincidence with the exact formulae 
developed recently by Cardy. 



1. INTRODUCTION 

It was a common belief until a very recent time 
that on two-dimensional (2D) lattices at perco- 
lation threshold pc there exists exactly one per- 
colation cluster New insight developed re- 
cently by Aizenman, who proved |^ that the num- 
ber of Incipient Spanning Clusters (ISC) in 2D 
critical percolation can be larger than one, and 
that the probability of at least k separate clus- 
ters is bounded A e"" < PL(fc) < e""' 
where a and a' are constants and L is a linear 
lattice size. We investigate by Monte-Carlo the 
number of spanning clusters in the critical bond 
percolation model on 2D square lattices. Using 
simple finite-size scaling of probabilities on a self- 
dual lattices of moderate size, we have determined 
with a good accuracy values of the probabilities 
P{k) = liniL^oo PL{k) for fc = 1, 2 and 3. 



twice the number of sites. Second, the self-duality 
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2. LATTICE 

We use in simulations a rectangular square lat- 
tice with linear size L, and exactly L sites and L 
bonds both horizontally and vertically. For clar- 
ity the example of lattice is shown in Fig. |l| to- 
gether with the dual lattice. The dual lattice has 
the same number of sites and bonds as the orig- 
inal lattice. This gives the possibility to keep in 
the finite lattices some properties of infinite lat- 
tice. First of all, the number of bonds is exactly 
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Figure 1. Example of lattice with L = 5. 



is valid for any finite lattice size. Third, the hor- 
izontal and vertical directions are equivalent. Fi- 
nally, one could expect trivial finite-size scaling 
for the bond percolation on such lattices. Indeed, 
probabilities are proportional to the inverse sys- 
tem "volume", i.e. cx l/L^ (Fig. | and Fig. |). 

3. ALGORITHM 

Each bond could be occupied with probabil- 
ity q = I — p and closed with the probability p. 
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Given the realization to each of 2 • bonds to be 
occupied formed the sample. Each sample could 
be exactly decomposed in clusters of connected 
occupied bonds. We are interested in the span- 
ning properties of such a clusters. Namely, what 
is the probability that a cluster connects the op- 
posite borders of square and what is the number 
of disjoint spanning clusters? It should be noted 



8, 12, 16, 20, 30, 32, 64 looking for aU events up to 
fc = 8. See for details the more full text [Q. On 
the Fig. H we plot the probability for simultane- 
ous occurence of more than one spanning clus- 
ter PL{k > 1) = J2k>i'^hk versus the inverse 
system volume A linear fit gives us the 

limiting value of P{k > 1) = 6.58 • lO'^ with 
error « 3 • 10^''. Fig. ^ shows the dependence 
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Figure 2. Probability of more than one Incipient 
Spanning Cluster. 



Figure 3. Probability of more than two Incipient 
Spanning Cluster. 



that the variance of probabilities is independent 
from the lattice size because probabilities are cal- 
culated as the expectation values of correspond- 
ing indicator functions. So, we should keep the 
number of samples independent of the lattice size 
and, therefore, the only parameter which controls 
the statistical errors is the number of samples M. 
Probabilities were calculated averaging M — 10^ 
samples and the error bars were defined from 100 
bins, each bin being the average over 10^ samples. 

4. NUMERICAL RESULTS 

We calculate the crossing probabilities tt^^ and 
7r„^ for square lattices with linear sizes L = 



of the probability of more than two clusters ver- 
sus 1/L^. The best linear fit gives the limiting 
value of P(fc > 2) = 1.48 • 10"'' with uncertainty 
2.1 • lO-'^. 

Actually, we observed in computations mostly 
the events of up to three simultaneous spanning 
clusters. We simulated at total about 10^° sam- 
ples of different sizes 8 < L < 64 and only one 
sample with 4 spanning clusters was detected. 
This event clearly not contradicts with our es- 
timate Q for P(4) « 10^^^. Single event is prob- 
able as one part in ten. 

A best linear fit to the logarithm of P{k) versus 
k^ gives a oc 1.61(7) and ^ oc 3 ± 1.5. 
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The crossing probability with periodic bound- 
ary conditions (PBC) is known to be larger and 
the value 0.63665(8) is computed in We esti- 
mate in computer simulations the probability of 
disjoint ISC clusters to be P(fc > 1) « 2.0(4) lO^^ 
and P{k > 2) « 1-4(5) 10~^ which are even 
smaller than for the case of free boundary con- 
ditions (FBC) considered by us. The finite size 
scaling of PBC is more complicated Q and fit of 
the results for only four lattice sizes (8,16,32,64) 
we simulated gives less accurate limiting values 
in comparison with the one for FBC. Surprisingly 
for us, the linear fit of logarithm of P{k) versus 
fc^ gives better accuracy then for the case of FBC: 
apBC « 1.915(1) and A w 4.26(3). The different 
behaviour of probabilities will be explained in the 
next section. 
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5. COMPARISON WITH EXACT RE- 
SULT 

Quite recently, John Cardy, using methods of 
conformal field theory, conjectured an exact form 
j5| of probability of k distinct Incipient Spanning 
Clusters in the limit of large aspect ratio r for 
rectangle (in our notations, FBC - free boundary 
conditions) and cylinder (PBC - periodic bound- 
ary conditions). It is possible to compare the data 
of our paper with the ratio, developed from for- 
mulae (2) and (3) of Cardy preprint 

logPFBcjk) ^ k{k-\) 
logPpBcik) fc2_i ' 

which is I = 0.8 for A: = 2 and f w 0.857 for 
k = 3. Our data gives 0.808(10) and 0.851(20), 
correspondingly. The discussion of that close 
agreement can be found in Cardy's preprint. The 
difference in the fc-dependence explains why our 
PBC data fits more closely to k^ than the FBC 
ones. Therefore, using the known /^-dependence 
of the probabilities, we can collapse all but fc = 1 
data on a single curve (Fig.^) versus a Cardy 
variable X which is XpBc = fc^ — 1/4 for pe- 
riodic boundaries and XpBC — k'^ — k/2 for free 
boundaries. The linear fit should give us an ap- 
proximation of the exponent aasymp in (2), which 
is 27r/3 = 2.094 by Cardy's conjecture. The linear 
fit to our data gives 1.91(20), which is, again, not 



Figure 4. Probabilities of fc = 1,2 and 3 simulta- 
neous clusters on lattices with periodic and free 
boundary conditions. Xpsc = fc^ — 1/4 and 
XpBC = fc^ — fc/2, correspondingly. 



so bad an approximation to Cardy's asymptotic 
value. 

It seems the fact that there are nonzero proba- 
bility (about half of percent) of the coexistence of 
infinite clusters of black (open) and white (closed) 
sites (bonds) in critical two-dimensional percola- 
tion. 
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